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(54) Real-time nitrogen oxides (NOx) estimation process 



(57) A process for providing a real-time estimation 
of in-cylinder NOx production rates during a combustion 
process in an internal combustion engine. The process 
employs a mode! that takes into account fundamental 
in-cylinder variables such as tennperature, pressure and 



oxygen content. The process can be used with calibrat- 
ed or uncalibrated internal combustion engines and can 
be used to capture the effects in variations of engine 
speed, engine torque. EGR, injection timing, and engine 
boost levels. 
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Description 
BACKGROUND 

5 [0001] This disclosure relates generally to diesel engines and more particularly, to systenns and methods for real- 
time estimation of NOx in exhaust gases generated by the diesel engine. 

[0002] Cun-ent production diesel engines generally do not employ any type of NOx estimation system. By omitting 
any meaningful real-time NOx estimation, it is possible that the diesel engine will perform sub-optimally due to the fact 
that other engine operational parameters cannot be dynamically adjusted due to potentially changing NOx levels that 

10 occur during the combustion process. 

[0003] Cun-ent production methods of estimating NOx, if employed, are generally based upon a lookup table that 
can only be generated for a fully calibrated engine, i.e.. once the exhaust gas recirculation (EGR), boost, fuel quantity, 
and timing controllers are designed. When using such techniques, real-time control of NOx is not possible because to 
control NOx, some of the actuators mentioned above are to be changed, thus changing the original values contained 

IS in the look up table. Moreover, the use of lookup table methods primarily involves steady state estimation of NOx. Thus, 
real-tlhie estimation is not provided by these methods and the use of these methods provides an inherent source of error. 

BRIEF SUMMARY 

20 [0004] Disclosed herein is a process for providing a real-time estimation of in-cylinder NOx production rates during 
a combustion process in an Internal combustion engine. The process comprises receiving a plurality of input values 
from a combustion chamber of a cylinder, wherein the Input values consist essentially of a mass amount of trapped air 
in the combustion chamber, a mass amount of tra|iped fuel In the combustion chamber, a total mass of species gen- 
erated within a combustion zone as a fraction of total in-cylinder mass, and a correlation value from a lookup table 

25 derived from an average NOx temperature of a combustion zone in the combustion chamber, and calculating an output 
signal from the pluralrty of input values, wherein the output signal produces a real-time estimate of the NOx production 
rate from the cylinder. . 
[0005] In accordance with another embodiment, a process for providing a real-time estimation of in-cylinder NOx 
production rates during a combustfon process in an internal combustion engine comprises calculating an average 

30 temperature value in a combustion zone of a cylinder combustion chamber con-elating the average temperature value 
to a correlation value in a look-up table; and estimating the in-cylinder NOx production rate from a cateulation including 
the con^elation value, an amount of trapped air in the combustion chamber, an amount of trapped unbumed fuel in the 
combustion chamber, and a total mass of species generated within a combustion zone as a f racoon of total in-cylinder 



35 [0006] The above described and other features are exemplified by the following figures and detailed description. 
BRIEF DESCRIPTION OF THE DRAWINGS 

[0007] Referring now to the f igures, which are exemplary embodiments and wherein the like elements are numbered 
alike' 

[0008] Figure 1 is a graphical illustration showing steady state and equilibrium cateulations for NOx mole fractions 
as a function of crank angle degrees; 

[0009] Figure 2 is a ptetorial illustration of combustion zone and an excess airzone vinth unbumed fuel In acombustion 
chamber cylinder at peak In-cylinder average temperature condittons; 
45 [0010] Figure 3 Is a block diagram of a NOx estimation model In a Matlab^SlmulInk fonnat; 

[0011] Figure 4 is a graphical Illustration of model predictions and measurement of NOx flow rate as a function of 
test number for a 1 .3 liter engine running at 1 .500 revolutions per minute (rpm); 

[0012] Rguie 5 is a graphical illustration of model predictions and measurement of NOx flow rate as a function of 
test number for a 1 .3 liter engine running at 1 ,800 rpm; 
so [0013] Figure 6 Is a graphical illustration of model predictions and measurement of NOx flow rate as a function of 
test number for a 1 .3 liter engine running at 2,400 rpm; 

[0014] Figure 7 is a graphical illustration of model predictions and measurement of NOx flow rate as a function of 
test number for a 1 .3 liter engine running at 2,800 fpm; 

[0015] Figure 8 is a graphical illustration of model predictions and measurement of NOx flow rate as a function of 
55 test number for a 1 .3 liter engine running at 3,000 rpm; and 

[0016] Figure 9 is graphical illustration showing error comparisons for the data illustrated in Figures 4-8. 
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DETAILED DESCRIPTION 

rnni-ii Disclosed herein is a process for providing real-time estimation of nitrogen oxide (NOx ) production rates 
Kacombl^^^^^^^^ 

As used heTSn the tem, "real-time" is defined as an instantaneous estimation of the amount of NOx currently being 
generSd au^^^^^^^^ of time within each cylinder. The process for providing rea«.ne NOx -t.ma.on empb^^^ 
model that takes into account fundamental in-cylinder variables such as temperature, pressure, oxygen content, and 
othJ ike i^r^e^^ the inK:ylinder model provides a robust method for NOx estimation that is surtably ^ph- 

^^nnd^^atefor^ienLranJofoperatingspeedsandl^^ 

jSm ratS^boost pressures ranging from no boost to the maximum boosting possible (below over-boosting), and 
iSTa vaTa^^ons whLsimultaneous^^^^^^ 

SvJ me s!mp^c^ of the in-cylinder model minimizes computational overhead on the power tram control module 
and can be used to eliminate the use of an expensive engine out NOx sensor. 

S)18J To define the In-cylinder model, engine combustion and emissions data were anah^ed to '^ent J/ t^e in 
Sender parameters having the greatest effect on NOx emisstons fom^aBons. Parameters were then identified on a 
S^^Cuon c^ctl^^^^^ Which altows the process to predict the transient behavtor and em-.^.ons cyde-to^cle van- 
ations Advantageously, the in^linder model can be easily Integrated with combustion models^ 
Sjl^l in order to develop the real-time emissions estimation system, some assumptions and sin^pllficatlons to the 
Son! lormt on process were made. For example, as shown in Figure 1 . a comparison of the NOx mole fraction 
fo^arnWs^.^ .0 the equilibrium NOX throu^ 

l3t„tr«tioi^^^^^^ lower than the steady state values after a peak temperature point. This obsenration is believed 
"Cra^Sl'othe n^^^^^^ kinetics for the NOx reactions. However. H was also observed that duri^ 

heLZsfonstrl^^^^ 

^etempcra^^^^^^^ 

lilhTzet "ox concentration is about constant at about 90% of the maximum fom,ed at the peak «-Pe^«^- 
roSoi Another assumption and simplrficatfon to the emissions fomiatlon process « shown in Figure 2. As show^ 
Lfinder has been SivWed into two zones: a combustion zone and an excess a.r zone The excess air zone « 
dSneTS the region of the cylinder with unbumed fuel and air. The combustion zone is defined by the Presence ^f 
the cher^fcal reactions and fSmation of NOx species as weD as other emissions. The combustion zone can be char- 
a^er^^Shrngthe bumedfuel at peaktemperature reacting 

an tSess amount of nitrogen due to the re-circulated exhaust (i.e.. EGR). As such, the oxygen percentage in the 
Shrs!Sl^rnLo)^ninthefreshair.Atpeaktemperatureconditions.i^^ 

teS^Sn threquilibrium Nbx fomiations with the steady state measurements. Therefore, the chem«al eq"-'^""^ 
SSonscrbe perfomiedforthe combustion zone atthe time of peak-m^^ 
tTpCie an estimate of the amount of NOx produced. The cateulations are <*°"^ 

ternoerature not the overall in-cylinder temperature. The chemical equilibnum cateulations Inc udemlnlmKing the total 
<^fre^ei;^Xallofthe2e^^ 

^^er^iS- refe^to gaseous components f omied in the mixture (combustton zone) including but not limrted 
to ^he^trSS oxi<^i lased upon this obsen^ation. a look up table can be made based on the <*emK«l equilibnin^ 
Sl^XSrr^latetotheLragetemperatureconditionforth^ 

f^ir?i?^Smbustlon zone is a non-unrfomi temperature distrtoution area, which can result in "0"-"""^°'^ f^'^" 

s^onsom.ation.lnthepresentmodel.itisassumedthattheaveragetemperatureofthecombus^^^^^^ 

NOX producing temperature. However, in most cases, this average NOx produang t^^^ff^^^ "^^"t from the 

acwal average temperature due to the non-linear behavior of emissions production as a function <>' temperattire 

rS ?3nTcr.elatlon of the average NOX producing temperature was dete^^^^^^^^ 

L easily implemented, yet can be accurate enough to capture the various effects. The estimatK^n model ctependant 
on the tapped airmass!^the cylinder pressure at maximum te-Perature. andthe^.nder pjessu^^^^^^^ 
Usina these input variables, the in-cylinder model output provides a real-time estimate of the NOx flow rate from the 
^Ser ^is m'^ mL only once^rcycle afterthe maximum temperature has occurred, whfch provides a significant 
advantage over chemical reaction integraUon models in real-time implementation. .„^^„ ^„ can b- 

[00231 The in^linder correlation model for providing real-time NOx estimation dunng combustion cycles can be 
determined from the following mathematical relationship as represented by fonnula (I). 
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Wherein is the trapped air mass (expressed in kiiograms), is the trapped mass at standard temperature and 
pressure with a volumetric efficiency of 1 (expressed in kilograms). is the cylinder pressure at maximum tem- 
perature (expressed in kilopascals), Ptdc is cylinder pressure at top dead center (expressed in kilopascals), and a 
(expressed in degrees Kelvin), b, and c are constants. It should be noted that T^^ represents the average NOx pro- 
ducing temperature in the combustion zone as previously described and shown in Figure 2 rather than a bulk temper- 
ature of the burning gas. In operatton. the Tb„„ value is then con-elated to a look-up table that includes static values 
based on the chemical equilibrium cateulations. Since the look-up table is based on chemical equilibnum calculations 
and as such, the values contained within the table are statte. The values in the look-up table do not change and are 
independent of actuator changes, engine types, and does not require a steady state estimation of NOx. It is further 
noted that constants a. b, and c did not show a significant change when increasing the number of test points used for 
idenirication. 

.. .uk <.>. ^^^^ ~«~.r>KIIKu nf <^ar>tiirlnn thA Hiff«>rAnt AffA^^S that can 

[0024] I ne mocei provioes excepuoiitii uwuiauy wiui mc wm»«»w...»^ w. ^^.-....9 

cause changes in NOx productton rales. Although the in-cylindermathematfcal model shown In Fomnula (1) is prefen^ed, 
a two-dimensional took up table can also be employed to provide similar results. The two-dimensional look-up table 
could use the mass ratios and pressure ratios as the input values for providing the conrelation. Validation of the model 
was conducted using a 1 .3 liter diesel engine with the following inputs recorded as shown in Table 1 . The NOx flow 
rate was estimated under a variety of conditions using the in-cylinder model and compared to actual measurements 
taken under the same conditions. Engine speed was varied to 1 500. 1 800, 2400, 2800 and 3000 revolutions per minute 
(RPM); toads (BMEP) were varied to 0. 200. 400. 600. and 800 kilopascals (kPa); EGR valve duty was vaned to 0, 30 
and 60 percent. Injection timing was varied at timed (calibrated), +5 degrees, and -5 degrees. 



Table 1. 



CONTROLLING 


RECORDED VARIABLES 


RPM (engine speed) 


Torque (engine torque) 


BMEP (brake mean effective speed) 


IMEPs (individual cylinders mean effective pressures) 


EGR DTY (EGR valve duty cycle) 


MAIR (air flow rate) 


SOI (start of Injection angle) 


MAP (intake manifold absolute pressure) 




IMT (intake manifoM temperature 




C02Jn, C02_out (intake and exhaust manifold carbon dioxide 
concentration) 




NOx_eng (exhaust NOx concentratfon (parts per million)) 




02_eng (exhaust oxygen concentration) 




A/F bench (emissions bench air to fuel ratio) 




SOI jpWoX (pilot start of Injection angle) 




fueLmaIn (main fuel injection desired volumetric flow rate) 




fuelj>ilot (pilot fuel injection desired volumetric flow rate) 




smoke (Bosch smoke number) 




EGR % (percentage of EGR to fresh air) 
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[0025] Figure 3 shows implementation of the in-cyllnder model in a Simulink function block format. The inputs for 
estimating NOx output include (1) trapped air. fuel and EGR masses for each cycle, (2) exhaust gas concentration, i. 
e.. a total mass of species generated within a combustion zone as a fraction of total in-cylinder mass determined from 
a ratio of the combustion zone mass to the in-cylinder mass. (3) fuel f ractbn burned up to maximum temperature time, 
and (4) a chemical equilibrium calculatfon value based upon the T^^^ con-elation, which is dependent upon the trapped 
air mass, the cylinder pressure at maximum temperature, and the cylinder pressure at top dead center. 
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rO0261 The results comparing the estimated NOx output values versus the actual measurements are graphcaHy 
Sown in Rqu^s 4 8 RgJre 9 displays the full test matrix used and presents the overall correla .on of the rnodel wrth 
tt^axZ^J^^rS^^F^gure 9 also shows the residuals and the ern,r percer,tage. The results deariy indicate the 
the artual measu^^^ The predicted NOx values closely overiapped the actual measurements under 

s :^stronSion^^^^^^ 

cor^oared to *e re^fts obtained under most conditions. While not wanting to be bound by theonr. a possible reason 
Sesf rI^2leV h^^^ enors is either a) a result of low measured values, wherein smaU absolute error values may 
nLsi^L havrresutted in high em>r percentages such as when observing the residuals for these points, or b) the 
c^Tustior^rS^^^^^^ smoke numbers, which suggests that the model m«ht not be as accurate 

10 at hiah eauivalence ratios compared to lean combustion conditions. ^ . ^^«,k..«. 

AHtegeously. the in-cylinder model for providing real-time estimates of NOx production durmg a combus- 
tion cycle offers at least the following advantages: 

1 . Iteal-time "transient" model, although it is based on steady state data but it Is a function of the in-cylinder cycle- 
15 by-cyde parameters; 

2. Runs once per cycle, which minimizes the computational expense; 

3. Easy inlegration with combustion models that are used to estimate the in-cylinder combustion characteristics; 

4. Reduced effective parameters that affect the estimation of NQx production from diesel combustion; 

5 Use of the model can potentially eliminate the use of engine-out NOx sensors that may eventually be required 
for production diesel engines in order to meet future emission standards; 

6. Calibration of the engine is not required to employ the model; and 

7. The look-up table is based on chemfeal equilibrium cateulations. and as such, the table values are static and 

do not change. 

r00281 While the disclosure has been described with reference to an exemplary embodiment it will be "nde|J«>J 
KseSS^ in the^that various 

wih^rdeoartnq from the scope of the disctosure. In addition, many modifications may be made to adapt a parhcular 
riLatlTmteriaTS t^eSchings of the disclosure without depaning from the essential scope thereof. Ther^o^ 
ft feSdeSt the disclosure nofbe limited to the partfcular embodiment disclosed as the best mode contemplated 
l^^^^o^S^ls^^^re. but that the disclosure will include all embodiments falling within the scope of the 
appended claims. 
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A process for providing a real-time estimation of In^linder NOx production rates during a combustion process in 
an internal combustion engine, the process comprising: 

receiving a plurality of input values from a combustion chamber of a cylinder, wherein the '"f"^^^'"^^ 
SeSly Of a mass amount of trapped air in the combustion chamber, a nnass amount ot t^PP^^^ 
combustion chamber, a total mass of species generated v^rithin a combustion zone as a fraction of total in 
Xer mai^^ a correlation value fSm a lookup table derived from an average NOx temperature of the 
combustion zone; and 

cateulating an output signal from the pluraHty of input values, wherein the output signal produces a real-time 
estimate of the NOx production rate from the cylinder. 

The process according to Claim 1 . wherein the average NOx temperature value is detemjned ^^^PP^^ 
mass value, a cylinder pressure at top dead center, and a cylinder pressure at maximum temperature. 

The process according to Claim 1 . wherein the average NOx temperature value is calculated in accordance with 
the following mathematical relationship: 
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wherein Th„„ is the average NOx temperature, ma,, is a trapped air mass value, m^p is a trapped mass value at 
standard temperature and pressure with a volumetric efficiency of 1. is a cylinder pressure at maximum 
temperature, Pjdc » a cylinder pressure at top dead center, and a. b. and c are constants. 

4. The process according to Claim 1. wherein the average NOx temperature is determined at a pealc in-cylinder 
average temperature condition. 

5. The process according to Claim 1 . wherein the combustion zone comprises reacting burned fuel at a peak tem- 
perature with a stoichiometric amount of air and a rediculated exhaust gas. 

6. The process according to Claim 1 . wherein the correlaaon value comprises a chemical equilibrium calculation for 
NOx species produced In the combustion zone. 

» 7 The process according to Claim 6. wherein the chemical equilibrium calculation comprises determining a mole 
fraction of species produced in the combustion zone at the average temperature value by minimizing a total Gibbs 
free energy for ail species. 

8. A process for providing a real-time estimation of in-cylinder NOx production rates during a combustion process in 
25 an Intemal combustion engine comprising: 

calculating an average NOx producing temperature value in a combustion zone of a cylinder combustion cham- 
ber; 

correlating the average NOx producing temperature value to a correlation value in a look-up table; and 

estimating the in-cylinder NOx production rate from a cateulation Including the correlatton value, an amount 
of trapped air in the combustion chamber, an amount of trapped unbumed fuel in the combustion chamber, 
and a total mass of species generated within a combustion zone as a fraction of total in-cylinder mass. 

9. The process according to Claim 8. wherein the average temperature value is determined from a mathematical 
relationship comprising: 
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45 Wherein Tk^™ is the average temperature value, mair Is a trapped air mass value, m^p Is a trapped mass value at 

standard temperature and pressure with a volumetric efficiency of 1, Py^ax a cy«nder pressure at maximum 
temperature, Ptoc "S a cylinder pressure at top dead center, and a. b, and c are constants. 

10. The process according to Claim 8. is dependent on in from a mass air flow value, a cylinder pressure at top dead 
so center, and a cylinder pressure at maximum temperature wittiin a combustion zone. 

1 1 . The process according to Claim 8. wherein the correlation value comprises a chemfcal equilibrium calculation for 
the NOx species. 

55 12. The process according to Claim 1 0. wherein the chemical equilibrium calculation comprises determining a mole 
fraction of species produced in the combustion zone at the average temperature value by minimizing a total GibDs 
free energy for all species. 
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